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ABSTRACT 

We report the results of a spectroscopic search for debris disks surrounding 41 
nearby solar type stars, including 8 planet-bearing stars, using the Spitzer Space 
Telescope. With accurate relative photometry using the Infrared Spectrometer 
(IRS) between 7-34 /mi we are able to look for excesses as small as ~2% of 
photospheric levels with particular sensitivity to weak spectral features. For 
stars with no excess, the 3a upper limit in a band at 30-34 jum corresponds to ~ 
75 times the brightness of our zodiacal dust cloud. Comparable limits at 8.5-13 
correspond to ~ 1,400 times the brightness of our zodiacal dust cloud. These 
limits correspond to material located within the <1 to ~5 AU region that, in our 
solar system, originates from debris associated with the asteroid belt. We find 
excess emission longward of ~25 /im from five stars of which four also show excess 
emission at 70 /im. This emitting dust must be located around 5-10 AU. One 
star has 70 micron emission but no IRS excess. In this case, the emitting region 
must begin outside 10 AU; this star has a known radial velocity planet. Only 
two stars of the five show emission shortward of 25 /im where spectral features 
reveal the presence of a population of small, hot dust grains emitting in the 7-20 
/im band. The data presented here strengthen the results of previous studies to 
show that excesses at 25 /an and shorter are rare: only 1 star out of 40 stars 
older than 1 Gyr or ~ 2.5% shows an excess. Asteroid belts 10-30 times more 
massive than our own appear are rare among mature, solar-type stars. 
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Subject headings: (stars:) circumstellar matter, debris disks, Kuiper Belt, Aster- 
oid Belt, comets, infrared, extra-solar planets 

1. Introduction 

The infrared excess associated with main sequence stars (Aumann et al. 1984) is pro- 
duced by dust orbiting at distances ranging from less than 1 to more than 100 AU. In 
discussing these debris disks we often use the analogy of our Solar System, where dust 
is generated by collisions between larger bodies in the asteroid and Kuiper belts, as well 
as outgassing comets. As in the Solar System, we distinguish between: a) hot (200-500 
K) debris located relatively close to the star (<1 AU out to ~ 5 AU) and arising in 
a region dominated by refractory material, e.g. debris associated with an asteroid belt 
(Dermott et al. 2002); and b) cooler debris (< 30 - 200 K) located beyond 5 AU that 
may be associated with more volatile, cometary material, i.e. the analog of our Kuiper 
Belt. Observed dust luminosities range from L dust /L+ ~ 10~ 5 to greater than 10 -3 for 
main-sequence stars, compared to L& ast / L± ~ 10~ 7 — 10~ 6 inferred for our own Kuiper belt 
(Stern 1996; Backman, Dasgupta, & Stencel 1995) and ~ 10 -7 measured for our asteroid 
belt (Backman & Paresce 1993; Dermott et al. 2002). 

The FGK survey of solar-type stars is a Spitzer GTO program designed to search for 
excesses around a broad sample of 148 nearby, F5-K5 main-sequence field stars, sampling 
wavelengths from 7-34 /zm with the Infrared Spectrograph (IRS; Houck et al. 2004) and 24 
and 70 /zm with the Multiband Imaging Photometer for Spitzer (MIPS; Rieke et al. 2004). 
In the context of our paradigm, IRS is well suited to searching for and characterizing hot 
dust in an asteroid belt analogue while MIPS is better suited to studying cooler material in 
a Kuiper Belt analogue. The sample consists of two components — stars with planets and 
those without. The MIPS observations of planet-bearing stars are discussed in Beichman et 
al. (2005a) and of non planet-bearing stars in Bryden et al. (2006). This paper presents 
the results of the IRS survey of a sub-sample of 39 stars, 7 with planets, and the remainder 
without planets based on current knowledge from radial velocity programs. In addition, the 
paper includes IRS follow-up observations of 2 stars (one of which is planet-bearing) that 
were not originally scheduled for spectroscopy, but which were found to have a 70 /zm excess. 
Results for one star with a notable IRS spectrum, HD 69830, have been presented elsewhere 
(Beichman et al. 2005b) but will be summarized below in the context of the complete survey. 

The IRS portion of the program was designed with four goals in mind: 1) to extend the 
wavelength coverage beyond that possible with MIPS to look for dust emitting either short- 
ward or longward of MIPS's 24yum band; 2) to take advantage of IRS's broad wavelength 
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coverage to reduce systematic errors and thus, possibly, identify smaller excesses than previ- 
ously possible; 3) to constrain the spatial distribution of dust grains by assessing the range 
of temperatures present in a given debris disk; and 4) to search for mineralogical features 
in the spectra of any excesses. In this paper we describe the sample of stars (§2); present 
our calibration procedure and give the reduced spectra for each star (§3 and Appendix-1); 
discuss detections of or limits to IR excesses (§4); present models of the excesses in the 4 
stars with IRS and MIPS 70 um excesses and discuss the nature of their debris disks; discuss 
overall results and conclusions (§5 and §6). 



2. The Sample 

The details of the sample are presented in Bryden et al. (2006). In brief, the sur- 
vey consists of solar-type stars with main sequence spectral types ranging from F5 to K5, 
most of which are within roughly 25 parsecs. Some stars at larger distances were included 
because they were known to possess planets via radial velocity studies (Butler et al. 1999; 
Marcy et al. 2002). A sub-sample of 39 stars covering the full range of spectral types 
was selected for IRS observations. Table 1 describes 41 stars including follow-up observa- 
tions of two stars earlier found by MIPS to have excesses (Table 4) (Beichman et al. 2005a; 
Bryden et al. 2006). Of the 39 original survey stars, seven are known to have planets from 
radial velocity studies with one additional planet-bearing star (HD 128311) selected for 
follow-up due to a possible MIPS excess. We also report MIPS results obtained after the 
observations reported in (Beichman et al. 2005a; Bryden et al. 2006) for 9 stars in the spec- 
troscopic sample (Table 3); three of the stars observed with IRS have not yet been observed 
with MIPS due to the as yet incomplete scheduling of the GTO program. Of the 9 new stars 
observed with MIPS, none has a 24 or 70 /im excess. 



3. Observations and Data Reduction 
3.1. IRS Spectra 

We observed each star with the three longest wavelength modules of the IRS: ShortLo-1 
(SL1, 7.5-14.0 /im), LongLo-2 (LL2, 14.0-20.5 fim) and LongLo-1 (LL1, 20-34 /im). The 
basic observing sequence and associated data reduction have been described in Beichman et 
al. (2005b) and are reviewed in the Appendix. In summary, we have used the fact that the 
majority of the sample shows no excess at MIPS wavelengths (24 or 70 /im) or in an initial 
examination of the IRS data to derive a "super- flat" to improve the relative and absolute 
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calibration of all the IRS spectra. After calibration with respect to photospheric models 
fitted to shorter wavelength (A < 5/xm) photometry, the agreement between the Spitzer/IRS 
data with IRAS and Spitzer/MIPS photometry is excellent: IRS/IRAS(12//m) = 0.96±0.01; 
IRS/IRAS(25^m)=0.93±0.01; and IRS/MIPS (24/im)= 0.99±0.03. 

To look for IR excesses we first examined differences between the flattened, calibrated 
spectra and photospheric models appropriate to each star (see Appendix-1 and Bryden et 
al. (2006)). For convenience, we define two "photometric bands" useful for isolating ei- 
ther the silicate features (8.5-13 /jm) or a long wavelength excess (30-34 //m). For all 
41 stars in the sample, the average level of the deviation from an extrapolation of the 
photospheric emission, (F u (Observed) — F u (Photosphere)) / F u (Photo sphere) , has values of 
-0.007 ± 0.067 in the 8.5-13 fim band and -0.006 ± 0.066 in the 30-34 jum band. Only 
HD69830 (Beichman et al. 2005b) stands out as having a significant excess from this first 
look at the data. 

However, a closer examination of the individual spectra showed that in most cases the 
data for an individual IRS module differed from the extrapolated photospheric emission in a 
systematic way that could be attributed to a simple scale or offset error. Without worrying 
about the source of the scale error (due to either the photospheric extrapolation or residual 
calibration errors in the IRS spectra), we derived a new calibration factor for each star and 
for each IRS module using the first 10 points of each module, effectively pinning the IRS 
spectrum to the photosphere at the short wavelength end of each of the three IRS modules: 
SL1 from 7.5 to 8.0 jum; LL2 from 14.0 to 19.8 /im; and LL1 from 19.5 to 21 fim. Values 
of this calibration factor deviated from unity on a star- by-star basis by less than 10%. This 
technique produced much smaller residuals showing no significant deviations from zero over 
the entire IRS wavelength range (Figures-la and b) for the vast majority of the sample. 
The dispersion in the deviation from a smooth photosphere was reduced from ~ 7% to 1.3% 
(8.5-13 /im) and 1.7% (30-34 jttm). 1 A longer wavelength excess beyond the first ten data 
points was not affected by the choice of the number of points used to set the normalization. 
Using this technique, we found no deviation from a stellar photosphere for the majority of 
stars. However, as described below, we did find clear evidence of an excess rising up at 
wavelengths longward of ~ 25 /im in five stars and a silicate feature appearing longward of 
8 jum in two stars. 

Unlike the other stars in the survey, HD101501 showed anomalous structure at the few 



: We also examined a star-by-star normalization factor derived using the entire module, not just the first 
10 points. This choice resulted in a similar reduction in dispersion, but had deleterious effects on the few 
sources identified with excesses rising to longer wavelengths. 
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percent level in its spectrum which could not be removed using the normalizations described 
above. The module-wide deviations from flatness are suggestive of uncorrected flat-field 
errors, not astrophysical effects. We put an upper limit of 5% of the photospheric emission 
on any excess for this star. One or two stars, e.g. HD76151 (Figure 4), showed small residual 
offsets of a few mJy between the spectrum and the photosphere in the 8-14 /im band where 
the star is brightest. As these "features" are at a level less than the ~1.5% dispersion in the 
excesses averaged over the entire sample, we do not regard them to be significant. 

After flattening and normalizing the IRS spectra as described above, we estimated 
the fractional excess (F u (Observed) — F u (Photosphere)) / F v (Photo sphere) . Figure 1 gives 
histograms of the fractional excess estimated in 8.5-13 and 30-34 /im "pass bands". In 
assessing the significance of an excess we looked at the internal uncertainty in the flux density 
measurement of a given star and the fractional excess relative to the ~ 1.5% dispersion of 
the entire sample (Table 2). Typically, the uncertainty in the fractional excess relative to 
the population as a whole is more important in assessing the reality of a spectral feature. 
There are a number of stars which appear to have a significant excess when looking only at 
the internal uncertainties in flux density, but which are not so impressive when compared to 
the dispersion in the overall population. If we require any excess to be more significant than 
3 times the population averaged dispersion (1.3% 8-13 /iin and 1.7% 30-34 /im, Figure 1), 
then based on the data presented in Table 2 we can claim statistically significant 30-34 /im 
excesses for five stars (HD 7570, HD 69830, HD 72905, HD 76151, HD 206860). One star, 
HD 69830, shows a highly significant silicate feature starting around 8 /im. Four of these 
sources, but not HD 69830, also show a 70 /iin excess. In the case of HD 72905 which also 
shows a 70 /im excess, its 8-13 fim flux density excess is highly significant relative to internal 
uncertainties (12<r), its fractional excess (2.8%) is modestly significant (2.2<r) relative to 
the dispersion in the population, and the spectrum (discussed below) closely matches that 
expected for silicate emission. Thus, we regard the 8-13 /im excess of HD 72905 as real. 

Figure 2a-f shows the IRS spectra for six sources (HD 128311, HD 216437, HD 7570, HD 
72905, HD 76151, HD 206860) plotted in three different ways: a) their calibrated spectrum; 
b) their excess relative to photosphere after normalization with respect to the first 10 points 
of each module, and c) the fractional amount of the excess relative to the photosphere after 
normalization. The dashed lines in the bottom panel show an estimate of the 2a dispersion 
in the deviations from the photospheric models. Fractional deviations between these lines 
should be regarded with skepticism. The first two stars (HD 128311, HD 216437) show no 
apparent excess and are representative of the bulk of the sample. The last four stars (HD 
7570, HD 72905, HD 76151, HD 206860, plus HD69830 shown in Beichman et al. 2005b) 
have excesses based on statistical criteria defined above. 
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3.2. New MIPS Data 

MIPS observations for nine stars not previously reported 
(Beichman et al. 2005a; Bryden et al. 2006) are summarized in Table 3. The 24 and 70 
/iin photometry were derived using processing steps detailed in Bryden et al. (2006) and is 
based on the DAT software developed by the MIPS instrument team (Gordon et al. 2005). 
For consistency, we use the same analysis tools and calibration numbers as were adopted by 
Beichman et al. (2005a) and Bryden et al (2006). 

At 24 /jm, we carried out aperture photometry on reduced images using an aperture of 6 
camera pixels (1 pixel=2.5" at 24 /im) in radius, a background annulus from 12 to 17 pixels, 
and an aperture correction of 1.15. The flux level is calibrated at 1.047 /iJy/arcsec 2 / (DN/s), 
with a default color correction of 1.00 appropriate for sources warmer than 4000 K at a 
weighted-average wavelength of 23.68 /im. 

At 70 /im we used images processed beyond the standard DAT software in order to 
correct for time- dependent transients, corrections which can significantly improve the sensi- 
tivity of the measurements (Gordon et al. 2005). Because the accuracy of the 70 /im data 
is limited by background noise, we used a relatively small aperture of just 1.5 pixels (1.5x9") 
in radius and a 4 to 8 pixel radius sky annulus. This aperture size requires a relatively high 
aperture correction of 1.79. The flux level is calibrated at 15,800 /iJy/arcsec 2 /MIPS_70_unit, 
again with a default color correction of 1.00. 

The 24 and 70 /im data are presented in Table 3 which also includes a comparison with 
the predicted level of photospheric emission based on models fitted to 2MASS and shorter 
wavelength data (Bryden et al. 2006). None of the new objects shows an excess at 24 /im. 
However, since accurate 2MASS data are not available for most of these bright (K s < 6 
mag) stars due to saturation effects, a significant part of the dispersion for the stars studied 
herein is probably due to the uncertainties in the photospheric extrapolation. We define the 
parameter %70 = (F u (Observed) — F u (Photo sphere)) / a at 70 /im where a is a combination 
of internal and calibration uncertainties as discussed in detail in Bryden et al. (2006). Using 
%70 as a criterion, none of the new objects show a significant 70 /tm excess. 

The 30-34 /tm excess toward the five stars described above (HD 7570, HD 69830, HD 
72905, HD 76151, HD 206860) have been carefully compared with MIPS observations at 
24 /im (Bryden et al. 2006). A sample of 87 stars in the FGK survey selected for hav- 
ing no 70 /xm excess (MIPS-70/Photosphere<1.5) has a 1-a clipped average of the MIPS- 
24/Photosphere ratio=0.978±0.001 with a dispersion of 0.056 (Bryden et al. 2006). HD 
69830 is clearly detected as a strong MIPS-24 excess with a MIPS-24/Photosphere ~ 1.5. 
HD 76151 and HD 206860 have MIPS-24/Photosphere ratios of 1.00 and 0.96 and show no 
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evidence for 24 /im excess. HD 72905 and HD 7570 represent intermediate cases, however, 
with MIPS-24/Photosphere ratios of 1.064 and 1.053. While the nominal dispersion for the 
24 jum data would suggest that these MIPS-24 excesses are not significant (< 2a), further 
analysis of the MIPS data (Kate Su, 2005, private communication) suggests that the com- 
bined photospheric and MIPS uncertainties may be closer to 2% in which case these modest 
excesses may turn out to be real. In a similar vein, Bryden et al. (2006) have identified a 
weak general trend (1.5 ~ 2a) for stars with 70 /im excesses to have slightly higher MIPS- 
24/Photosphere ratios than for stars without a longer wavelength excesses. If the weak MIPS 
excesses hinted at here can be made more significant with careful calibration and validation 
using IRS, then it may be possible to look for weak, warm excesses in a much larger sample 
of objects than would be possible with IRS alone. 



4. Results 

4.1. Statistics of detections 

We detected significant IRS excess emission in five sources including HD 69830 (Be- 
ichman et al. 2005b). Of these five stars, all except HD69830 also have a MIPS 70 /im 
excess. Two stars, HD 128311 and HD 206860, were included explicitly for IRS follow-up 
observations because they had MIPS excess emission. HD 128311 shows no IRS excess 2 
while HD 206860 does show an excess at 30-34 /im. Out of the 39 stars in the initial sur- 
vey, we have identified 30 jum IRS excesses around four stars (HD69830, HD72905, HD7570, 
HD76151), or 11±5%, which is consistent with the fraction of stars with excesses found at 
70 /im (Bryden et al. 2006). Including the two follow-up stars (HD206860 and HD 128311) 
brings the detection rate to 5 of 41 stars or 12±5% at 30/iin. None of the eight planet-bearing 
stars have 30 /im excesses, a smaller fraction than seen with MIPS 70 /im excesses, < 12% 
vs. 24±10% (Beichman et al. 2005a). The result is suggestive of an absence of warm dust at 
levels detectable by IRS, but is not inconsistent with the expected fraction given the small 
number statistics. 

If we apply a similar analysis to the 8-14 /im portion of the spectrum, the corresponding 
statistics give detections of just two objects (HD69830 and HD 72905) out of the 39 stars 
(with and without planets) in the initial survey, or 5.1±3.6%, and zero of eight planet-bearing 
stars. Removing HD 72905 from the overall statistics because it is quite young, 300 Myr 



2 The 70 /zm detection of HD 128311 (Beichman et al. 2005a) is marginal and awaits confirmation in a 
pending MIPS observation. 
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(Spangler et al. 2001) based on its membership in the Ursa Major moving group, reduces 
the statistics to just 1 short wavelength excess out of 38 mature stars (> 1 Gyr as discussed 
in Bryden et al 2006) or <2.5% in the initial survey. These data confirm the rarity of short 
wavelength excesses compared with ones at longer wavelengths. This result was first noted 
on the basis of IRAS and ISO data (Aumann & Probst 1991; Mannings & Barlow 1998; 
Fajardo-Acosta, Beichman & Cutri 2000; Laureijs et al. 2002) and reconfirmed with MIPS 
data (Beichman et al. 2005a; Bryden et al. 2006), but is extended here to lower levels of 
emission, i.e. 2% relative excess for IRS vs. 6-10% for IRAS and ISO. In the next section we 
convert these observational limits into a limit on the fractional luminosity of dust in these 
systems, L dust /L*. 



4.2. Limits on the Fractional Disk Luminosity, L dust /L+ 

A useful metric for the limits on dust surrounding these stars is which is related 

to the fractional flux limit of an excess relative to the Rayleigh Jeans tail of the star's 
photosphere (Bryden et al. 2006): 

-^dust -^dust hvTl ^ 



F-k ~ U fcT d 4 ust (eW^ st _i) 
From which it is easy to show that 



Ld F dast e Xd — 1 (T d , 



L* Fjf Xd 

where F dust = F v (Observed) — F u (photosphere) . At the peak of the black body curve 
x d = hu/kTd has a constant value of 3.91, corresponding to T d = 367 K at 10 fim. At this 

wavelength L d jU = 0.00345 (^^j F d /F*. In Table 2 and Figure 3, we evaluate L dust /L* 
for each star using the appropriate effective temperature and luminosity (from SIMBAD, 
Table 1) and its measured fractional excess, F dns t/F*, or in the case of an upper limit, 3a pop 
where a pop is the dispersion in fractional excess averaged over the whole sample. We calculate 
these dust limits in the 8.5-13 /im and the 30-34 /xm portions of the IRS spectrum assuming 
Td = 367K and T d = 115K at 10 and 32 /im, respectively. This definition Ld us t/£* assumes 
that the emitting material is all located at the location where the peak of the T d blackbody 
matches the wavelength of observation. More emitting dust, and higher values of L dust /L+, 
would be required for material located substantially interior or exterior to this point. 
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The 3a limits on Ldust/At at 8.5-13 /im and 30-34 /im have 2a clipped average values of 
L dust /L+ = 14 ± 6 x 10~ 5 and 0.74 ± 0.31 x 10" 5 respectively (Table 2). In comparison with 
our solar system, L dust /L^ ~ 10~ 7 (Backman & Paresce 1993; Dermott et al. 2002), the IRS 
results set limits (3a) on warm (360 K) dust peaking at 10 of ~ 1,400 times the level of 
emission in our solar system. For cooler dust (~ 120 K) peaking at 30-34 /im, the 3a limit 
corresponds to ~74 times the nominal Ld nst / L± of our zodiacal cloud. 

For objects with detections in the IRS bands, we can determine L dust /L^ explicitly by 
integrating over the data between 10-34 and extrapolating using the models discussed below 
out to 70 /im. These values are given in Table 4 and are compared with the 70 /zm-only 
estimates obtained using equation (1). Including the IRS spectra increases the estimated 
Ldust/L* by about a factor of two relative to the 70 /im-only estimate. For the nine stars 
with new upper limits to an excess at 70 /im (Table 3), we also report limits on L dust /L*. 

5. Discussion 

5.1. Models of Disks around Solar Type Stars 

We fitted the IRS spectra and MIPS 70 /im photometry using a simple model of opti- 
cally thin dust within one or two flat dust annuli centered around the star. As a guide to 
subsequent modeling, we first modeled the emission using blackbody grains at a single tem- 
perature in equilibrium with the parent star. From the ratio between the 30-34 /im and the 
70 /xm flux densities for each star we found (Table 5) grain temperatures in the range 75-100 
K with equilibrium distances of 8-14 AU for the four stars. Even though the uncertainties 
in the flux density ratio are quite large (~ 25%), the temperature and location of the grains 
are reasonably well defined, ±1 AU and ±15 K, given the assumption of blackbody grains. 

We calculated the power-law temperature profiles, T(r) = T o (L/L ) a (r/r o ) /3 , for these 
grains as a function of stellar luminosity and distance from the central star by solving 
for the grain temperature at which the integrals over the absorbed stellar radiation and 
emitted thermal radiation were equal. The equilibrium temperature and power-law vari- 
ation of temperature with distance and stellar luminosity depend on the assumed grain 
size and composition (Backman & Paresce 1993). As described below, we use 10 and 0.25 
/im silicate grains (Draine & Lee 1984; Weingartner & Draine 2001) for which the we ob- 
tained the following numerical relationships for large and small grains, respectively: T(r) = 
255 K(L/L Q )°- 26 (r/AU)-° Ag and T(r) = 362 K(L/L Q )°- 29 (r/AU)-° A1 . These calculated co- 
efficients and power-law constants closely follow analytical results (Backman & Paresce 1993). 

We then calculated the dust excess by integrating over the surface brightness of a 



-10- 



disk between Ri and R 2 , with F U (X) = J r (\)(r /r )~ p B u (T(r))rdr . In calculating the 
wavelength-dependent dust optical depth, r(A), we adopted grain emissivities Q a bs appropri- 
ate to small (0.25 /im) or large (10 /im) grains (Draine & Lee 1984; Weingartner & Draine 2001). 
The integrated surface density in the wedge-shaped disk expected for grains dominated by 
Poynting-Robertson drag is roughly uniform with radius, p = 
(Burns et al. 1979; Buitrago & Mediavilla 1985; Backman 2004). We examined a num- 
ber of other cases with < p < 1 that would reflect different dust dynamics, but did not find 
results that were substantially different from those for p — 0. Results of the model fitting 
are shown in Figure 4 and Table 5 and are discussed below. 

Two other Spitzer studies have recently been published with similar results to those 
found here. Hines et al (2006) identified a much younger (~ 30 Myr) star, HD 120839, with 
an IRS spectrum similar to that of HD69830. Reporting on early results from the FEPS 
survey (Meyer et al. 2004), Kim et al (2005) found a number of stars with 70 /im excesses 
with, in some cases, an accompanying excess in the IRS bands. The model results concerning 
the amount of dust and its orbital location are similar to those presented here. 

5.1.1. Single Ring Models for HD 76151, HD7570, HD 206860, and HD128311 

For HD 7570, HD 206860 and HD 76151, which show no statistically significant, short 
wavelength excess, we used only a single population of large, 10 /im, amorphous, silicate 
grains (Draine & Lee 1984; Weingartner & Draine 2001), fitting emission from a single an- 
nulus to 89 data points longward of 21 /im (just longward of the last point used for the flux 
normalization of the LL-1 IRS module) and including 70 /im. By varying To, R± and R2 we 
were able to minimize the reduced x 2 to values between 0.5-0.8 with 77 — 3 = 74 degrees of 
freedom. 

The models for HD 76151 and HD 7570 require grains with a narrow spread of temper- 
atures around 100 K and located in the range of 7-10 AU. Assuming from 10 /im grains, the 
predicted 850 /im flux density from the HD 7570 disk is well below the 1 mJy limit from 
SCUBA (J. Greaves, private communication, 2005). The predicted 850 /im flux density from 
a disk containing larger (~ 100 /im) grains this system would be approximately 1 mJy, just 
at the SCUBA limit. We cannot distinguish between the two grain sizes using only Spitzer 
data. HD 206860 has cooler grains (75 K) than the other stars with material located around 
14 AU. 

Mass estimates are notoriously tricky to derive given uncertainties in grain sizes, but 
these three systems suggest dust masses around 1-2.6 xlO~ 6 M ffi (Table 5) for a silicate 
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grain density of 3.3 gm cm -3 . Extrapolating this estimate using the —3.5 index power law 
appropriate for a distribution of sizes from a collisional cascade (Dohnanyi 1969) up to a 
maximum size, R ma x, of 10 km yields total mass estimates of 0.024-0.080 >J R max /10km M e . 
Submillimeter observations of all these sources would further constrain the dust size and 
distribution and thus the total mass of the emitting material. 

Since the IRS data place only a limit on the emission shortward of 34 /zm for the planet- 
bearing star HD 128311, the single grain model can be used to show that the inner edge of 
the disk seen at 70 /xm must start beyond ~ 15 AU corresponding to material cooler than 
~ 50K. 



5.1.2. Two Ring Model for HD72905 

Since HD 72905 shows evidence for emission from small silicate grains starting around 
8 jum, we separately examined three models: 1) small grains over a broad range of dis- 
tances; 2) small grains in two annuli, one close to the star and one further away; 3) small 
grains in a close-in annulus and large grains in a more distant annulus. For the small 
grains we used 0.25 /zm grains composed of either amorphous silicates or crystalline forsterite 
(Malfait et al. 1998; Jager et al. 1998). Unlike the dramatic case of HD 69830, small crys- 
talline silicate grains provide a poor fit to the 8-14 /jm data for this star (Figure 4). Amor- 
phous silicate grains have a broad emission profile particularly on the short-wavelength side 
which appears in the HD 72905 spectrum but not in HD 69830 's. A model providing a good 
fit to the data incorporates amorphous silicate grains located in an annulus extending from 
the grain sublimation distance (0.03 AU corresponding to 1500 K) out to 0.43 AU where the 
characteristic temperature is 500 K. The spectrum suggests an upper limit of no more than 
10% of crystalline forsterite grains (by number). 

However, this small dust population cannot extend out continuously much beyond 0.4 
AU. Even with a density falloff proportional to r~ 0A to r _1 , the presence of any dust between 
1 and 10 AU produces far more radiation in the 20-70 /zm region than is observed in the 
combined IRS+MIPS data. Thus we reject model #1 above. Rather, the data longward of 
14 /iin suggest a second population of cooler material located beyond 10 AU. We examined 
models with either small grains (0.25 /im, model #2) or large grains (10/xm, model #3) 
emitting at temperatures less than 90 K. Both models #2 and #3 give reasonable fits to 
the observations with reduced x 2 ' s ~ 0.8. As Table 5 indicates, a much greater quantity 
of small grains emitting at cooler temperatures is required to fit the data compared to a 
model using larger grains. This difference is due to the greatly reduced emissivity of small 
0.25/im vs. large 10/im grains, e.g. at 70 /im, Q abs =0.008 vs. 1.84 (Draine & Lee 1984; 
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Weingartner & Draine 2001). 

A testable difference between the large and small grain models is the size of the emitting 
region. For small grains, the emitting annulus extends between 109 and 310 AU, correspond- 
ing to a diameter of 40" at the 14.1 pc distance to the star. By contrast, large grains cool 
more efficiently and are located closer to the star, < 16 AU, for a predicted diameter of only 
2". Since there is no evidence that the 70 /zm image of HD 72905 is extended, we favor the 
large grain model. The total amount of material in the HD 72905 disk is about a factor of 
two greater than in the other disks studied here. 



5.2. Nature of Disks 

5.2.1. Rarity of Emission from Warm Dust 

The vast majority of the stars studied here, planet-bearing or not, have no emission 
interior to ~ 10 AU at the levels to which IRS is sensitive, 75-1,400 times that of our 
zodiacal cloud. Since this result applies equally to stars with and without planets, there 
must be some general explanation for the absence of warm material around mature stars 
(> 1 Gyr). A long term decline in the amount of debris expected on theoretical grounds 
(Dominik & Decin 2003) and demonstrated with the envelope of the distribution of 24 jum 
excesses seen in A stars, implying an excess oc 15 ^ yr (Rieke et al. 2005), is consistent with 
the results of this survey. Modifying Eqns. (35)- (40) of Dominik and Decin (2003) for 
parameters appropriate to an asteroid belt located between 1-5 AU leads to the conclusion 
that after one billion years, La ust /L* in the IRS wavelength range should be < 10~ 6 and thus 
undetectable with IRS over a broad range of initial conditions. 

Three of the disks (HD 7570, HD 206860, HD76151) have an inner edge around ~ 10 
AU and may be only a few AU wide (Table 5). The cutoff at wavelengths shorter than 
~25 fim corresponds to temperatures < 75-100 K and location > 7 AU. Other stars ob- 
served in this or other Spitzer programs with 70 fxm excess that do NOT show an inner 
emission cutoff at 25 /iin seem to have some extenuating circumstance: Vega with evi- 
dence for a recent collision (Su et al. 2005); (3 Pic which is extremely young; Fomalhaut 
(Stapelfeldt et al. 2004), two young stars HD 72905 and HD 12039 (Hines et al. 2006) and 
HD 69830 (Beichman et al. 2005b) with a sharply bounded, interior dust belt originating, 
perhaps, in an active asteroid belt. If the observational trend for a general lack of material 
interior to ~10 AU continues as more Spitzer data accumulates, then it may be that this 
inner boundary to the "Kuiper Belt" is related to the ice-sublimation distance > 3L°' 5 AU 
where significant amounts of material could have accumulated in the protoplanetary disk 
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(Thommes, Duncan & Levison 2003). 

The amount of warm material around HD 72905 with L dust /L^ ~ 18 x 10~ 5 is well above 
values consistent with the simple models (Dominik & Decin 2003) even given the relative 
youth of this star (~ 0.3 Gyr). The excess may be attributable to a recent collisional event 
between planetesimals as has been suggested by Kenyon and Bromley (2004) and Rieke et 
al. (2005) for A stars observed at 24 /j,m. As discussed in Beichman et al. (2005b), the disk 
of hot, small grains around a mature star, HD 69830, stands out as exceptional in this survey 
and may similarly be the result of a rare transient event in an unusually massive asteroid 
belt. In the case of steady-state emission (rather than transient events), the brightness 
of an excess will be proportional to the first or second power of the number of colliding 
planetesimals, depending on whether the ultimate loss mechanism for the smallest grains is 
blowout by radiation pressure or PR-drag, respectively (Dominik & Decin 2003). Since PR- 
drag is likely to be more effective than blowout for the low luminosity stars considered here 
(Sheret, Dent & Wyatt 2004), a 74-1,400-fold brighter zodiacal cloud than our own would 
imply the existence of (or a limit to) an asteroid belt \^75 ~ VT400 ~ 8 — 35 times more 
massive than our own, depending on its location. The results of this IRS survey suggest that 
asteroid belts of such a large mass are rare around mature, solar type stars. 

As pointed out by an anonymous referee, there are alternatives to the transient, col- 
lision hypothesis for the origin of the debris disks around mature stars. Stirring up a qui- 
escent asteroid or Kuiper belt by planetary migration caused by either orbital interactions 
(Gomes et al. 2005; Strom et al. 2005) or to the effects of a passing star (Kenyon & Brom- 
ley 2004b) has been suggested as possible causes for analogues of the period of Late Heavy 
Bombardment in our solar system. A strong infrared excess could be a signpost of a similarly 
exciting period in other planetary systems. For the young stars HD 72905 (this survey) and 
HD 120839 (Hines et al. 2006) this analogy is quite plausible. On the other hand, for the 
older stars in this survey, including the >2 Gyr-old HD 69830, it is not clear that planetary 
migration is likely so long after the formation of these planetary systems. Confirmation of 
this hypothesis will require a more complete census of planetary systems than is at present 
available. 



5.2.2. Lack of Warm Dust and the Presence of Planets 

None of the planet-bearing stars show an IRS excess and only one, HD 128311, shows 
a MIPS excess. On the assumption that the gravitational influence of a planet at orbital 
semi-major axis, a p , and eccentricity, t p is given by r grav ~ yj e^a^ + 12R 2 Hill where the Hill 

Radius is defined as R H ai = a p (M P /(3M*)) 1/3 (Bryden et al. 2000) with M p and M* the 
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masses of planet and star. As shown in Figure 5, the region between Ruin = a — r grav and 
Rmox = a + r grav could be swept clean of dust due to sweeping up action by the planet 3 . In 
addition to more global effects leading to a general decline in dust content as discussed in 
the preceding section, gravitational effects of planets in relatively distant orbits (HD 39091, 
55 Cnc, 47 U Ma, HD 216437, and v And) might clear out the material interior to ~ 10 
AU in the absence of a secondary source such as a massive asteroid belt. Stars with close-in 
planets (r Boo and 51 Peg) may require another explanation (or an as yet undetected, more 
distant planet) for the lack of material in the 1-10 AU region. 



6. Conclusion 

We have used the IRS spectrometer on the Spitzer Space Telescope to look for weak 
excesses around main sequence, solar type stars. After careful calibration, the IRS spectra 
allow us to achieve relatively low values of L dust /L+ compared to previous studies. At 8-13 
/im the IRS data provide an average 3a limit of L dust / ' L± = 14 ± 6 x 10~ 5 or roughly 1,400 
times the level of our solar system zodiacal cloud for material in the 1 AU region of the target 
stars. Two stars, HD 72905 and HD 69830, have detectable emission at these wavelengths 
and in both cases there is evidence for emission by small dust grains. At 30-34 /im, we reach 
a more sensitive, limit on L dust /L^ of 0.74 ± 0.31 x 10 -5 , or roughly 74 times the level of 
our solar system zodiacal cloud for material in the 5 — 10 AU region which may represent 
the inner edge of the Kuiper Belt. Four stars (HD 7570, HD 72905, HD 76151, HD 206860) 
showed evidence for long wavelength IRS emission that appears to link smoothly to a 70 
/im MIPS excess. HD 69830 shows no evidence for 70 /im excess and thus for colder, more 
distant material. 

The lack of IRS excess emission toward planet-bearing stars may be due to the small 
size of the sample observed to date, to the effects of planets removing dust from the orbital 
locations probed at IRS wavelengths, or to more general processes that remove dust from 
regions interior to a few AU (Dominik & Decin 2003). We suggest that asteroid belts located 
between 1-10 AU and as massive as 8-35 times our own asteroid belt are rare around mature, 
solar type stars. 



3 Planet properties from Planet Encyclopedia at http://cfa-www.harvard.edu/planets 
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8. Appendix-1. Details of IRS Data reduction 

Standard IRS Staring mode observations were made with the Short Low Order 1 (SL1; 
7-14 /iin), Long Low Order 2 (LL2; 14-20 /iin), and Long Low Order 1 (LL1; 20-35 /im) 
modules. Each star was observed at two positions along the slit, called Nodi and Nod 2 . 
Data were processed by the Spitzer Science Center (SSC) to produce calibrated images of 
the spectrometer focal plane. The data presented here were processed with Version 11.0 of 
SSC pipeline (February 2005). The SSC SPICE software was used to extract spectra from 
the images. Because of the need for careful subtraction of the stellar continuum to detect 
a faint excess, additional steps were taken in producing the final spectra. These are now 
described briefly. 

The first few points at the beginning and end of the spectrum from each module were 
typically unreliable, as were a few bad pixels flagged in the SSC processing. These data 
were rejected. This effect was particularly noticeable at the long wavelength end of SL1 
(A > 14 /im). However, since the short wavelength end of LL2 overlaps SL1, there is no gap 
in the final spectrum. 

Longward of ~7 fim, the photospheres of solar type stars are smooth and do not differ 
greatly from a Rayleigh- Jeans blackbody curve (Kurucz 1992; Castelli & Kurucz 2003). We 
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used this fact to improve the pixel-to-pixel calibration in the extracted SSC/SPICE spectra. 
Of the 41 stars listed in Table 1, 28 stars showed no clear evidence for an excess at either 
24 or 70 /iin with MIPS or with the initial pass of IRS processing but did show consistent 
pixel-to-pixel deviations from a smooth photospheric model. 

For these stars we formed ratios of the extracted IRS spectra to the Kurucz model 
appropriate for the effective temperature and metallicity of each star. From the average 
of these ratios at each wavelength, we created a "super-flat" response curve for the SL1, 
LL2 and LL1 modules. Only those stars which did not show signs of an excess during 
the initial pass of the data were included in the super-flat. Shortward of 14 /xm,the SL1 
super-flat typically has values in the range 0.97-1.03, or, equivalent ly, deviations from unity 
response of ~3%; a few pixels have values deviating from unity by 10%. The values at each 
pixel of the SL1 super-flat have a dispersion around the pixel response (averaged over 28 
stars showing no signs of excess in IRS or MIPS data) of a pop ~ 2% and a mean ~ 0.4%. 
In the LL1 and LL2 modules, the typical super-flat values are in the range 0.98-1.02 with 
dispersions of a pop ~ 1 — 2% and a mean ~ 0.2 — 0.3%. We divided the extracted spectra by 
the appropriate "super-flat" to remove any residual calibration variations. The dispersion 
in the super-flat ( o pop) is clll indication of the limiting systematic noise in removing the 
photospheric contribution to the signal from these stars. 

Although the spectrum of HD 101501 shows no excess at the ~10% level in either the 
MIPS or IRS data, its spectrum did not behave in a consistent manner when the various 
normalization steps were applied to look for fainter excesses. Differential pointing issues 
or other time variable effects could have produced a spectrum deviating from the mean 
for non-astrophysical reasons. The module-wide deviations from flatness are suggestive of 
uncorrected flat-field errors, not astrophysical effects. As a result, this star was excluded 
from the generation of the super-flat. 

The comparison between the SSC spectra and the photospheric models fitted to shorter 
wavelength photometry (Bryden et al. 2006) shows absolute calibration gain factors of 
F^SSC) /.^(Photospheric Model)= 0.84 ± 0.01 (5% dispersion from star to star) in the 
LL-1 and LL-2 modules and 1.12 ± 0.014 (7% dispersion from star to star) in the SL-1 mod- 
ule. As discussed in §3.1 the agreement between IRS data and IRAS or MIPS observations 
is excellent with deviations from unity of 4±1% at 12 /im, 7±1% at 24 /im (consistent with 
a known calibration issue with IRAS data at this wavelength, e.g. Cohen et al 1992), and 
1±1% at MIPS-24 /im. 
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Table 1. An IRS Survey of F, G, and K Stars 



Star 1 


Name 


Sp. Type 


V(mag) 2 


Dist(pc) 


AOR 


MIPS Excess 


HD 693 


GJ10 


F5V 


3.51 


18.9 


4008448 


No. New 


HD 3795 


GJ799 


G3V 


3.86 


28.6 


4009216 


N/A 


HD 4628 




K2V 


6.41 


7.4 


4009984 


No 


HD 7570 


GJ55 


F8V 


4.07 


15.1 


4010240 


70 /mi 


HD 9826 5 


v And 


F8V 


4.09 


1.35 


4010496 


No. New 


HD 10800 


GJ67.1 


G2V 


3.71 


27.1 


4010752 


No 


HD 39091 5 


GJ9189 


G1V 


5.67 


18.2 


4015616 


No 


HD 43834 


GJ231 


G6V 


5.05 


10.2 


4015872 


No 


HD 55575 


GJ1095 


GOV 


4.42 


16.9 


4016128 


No 


HD 58855 


GJ9234 


F6V 


3.87 


19.9 


4016384 


No 


HD 69830 


GJ302 


K0V 


5.45 


12.6 


4016640 


24 fim 


HD 72905 


GJ311 


G1.5V 


4.93 


13.8 


4016896 


70 (im 


HD 75732 5 


p (55) Cnc A 


G8V 


5.95 


12.5 


4017152 


No 


HD 76151 




G3V 


5.74 


11.3 


4017408 


70 /im 


HD 84737 


GJ368 


G0.5Va 


3.77 


18.4 


4017664 


No 


HD 86728 


GJ376 


G3Va 


4.53 


14.9 


4017920 


No. New 


HD 95128 5 


47 UMa 


G1V 


5.1 


14.1 


4018432 


No 


HD 101501 


GJ434 


G8Ve 


5.43 


9.5 


4018944 


No 


HD 120136 5 


t Boo 


F6IV 


4.5 


15.6 


4021760 


No 


HD 133002 


GJ3876 


F9V 


2.46 


43.3 


4022272 


No 


HD 136064 


GJ580.2 


F8V 


3.11 


25.3 


4022784 


No 


HD 142373 


GJ602 


F9V 


3.61 


15.8 


4023040 


No 


HD 146233 




G1V 


4.56 


15.4 


4081920 


No 


HD 154088 


GJ652 


G8IV-V 


5.3 


18.1 


4029952 


No. New 


HD 166620 




K2V 


6.14 


11.1 


4023808 


No 


HD 168151 




F5V 


3.18 


23.5 


4024064 


No 


HD 173667 


GJ9635 


F6V 


2.79 


19.1 


4024320 


No 


HD 181321 


GJ755 


G5V 


4.89 


20.9 


4030208 


No 


HD 185144 




K0V 


5.93 


5.6 


4024576 


No 


HD 188376 




G5V 


2.82 


23.8 


4025088 


No 


HD 191408 


GJ783A 


K3V 


6.41 


6 


4025600 


No. New 


HD 196378 


GJ794.2 


F8V 


3.2 


24.2 


4025856 


No. New 


HD 203608 


7 Pav 


F8V 


4.4 


9.2 


4026368 


No 


HD 212330 


GJ857 


G3IV 


3.76 


20.5 


4027136 


No. New 


HD 216437 5 


p Ind 


G2.5IV 


6.06 


26.5 


4084224 


No 


HD 217014 5 


51 Peg 


G2.5IV 


5.49 


15.4 


4027648 


No. New 


HD 217813 




G5 


4.72 


24.3 


4030464 


No. New 


HD 222368 




F7V 


3.43 


13.8 


4028416 


No. New 


HD 225239 


GJ3002 


G2V 


3.28 


36.8 


4028672 


N/A 


Follow-up 


Observations 












HD 128311 5 


GJ 3860 


K0 


7.51 


16.5 


4083712 


70 (im 


HD 206860 


GJ836.7 


GOV 


4.62 


18.4 


12719872 


70 /^m 



IRS Excess 4 



Long 



Short, Long 
Short, Long 

Long 



Long 



'Two other stars, HD 190248 (AOR 4025344) and HD 209100 (AOR 4026625), were observed as part of this program using 
the IRS High Resolution mode. These stars will be discussed elsewhere. HD 33262 was not observed properly due to a pointing 
error and is not considered further. 

2 Spectral Types from SIMBAD. Visual magnitudes and distances as quoted in SIMBAD, typically from the Hipparcos satellite 
(1997). 
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3 Star either does or does not have 24 or 70 (im excess as observed by MIPS. 'New' means MIPS data were obtained 
subsequently to data presented in Beichman et al 2005a or Bryden et al 2006 and are described in Table 3. 'N/A' means star 
has not yet been observed with MIPS due to scheduling constraints. 

4 Star has excess as observed by IRS as described in this paper at either Short (8-13 /im) or Long (30-34 /im) wavelengths. 

5 Star has at least one radial velocity planet. 
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Table 2. IRS Excess in 8.5-13 and 30-34 /iin Bands 





Excess (mjy) 


Frac. 


Excess / 


Ldust/M 1 ) 


Excess (mjy) 


Frac. 


Excess/ 


Ldust/M 1 ) 


Star 


(8.5-13 ^m) 


Excess 


Dispersion 


(io- 5 ) 


(30-34 fim) 


Excess 


Dispersion 


(IO" 5 ) 


HD 693 


OA 0_L1 Q 

-2U.o±4.o 


-0.015 


-1.2 


<8.2 


0.3±0.8 


0.002 


0.1 


<0.4 


tiu o7yo 


-5.5±0.5 


-0.008 


n t 
-U. 7 


<13 


-1.9±0.6 


-0.026 


-1.3 


<0.7 


44 4J 4bZo 


i q o _i_n n 
4o.o±U.y 


0.009 


n T 
U. 7 


<26.8 


-0.8±1.3 


-0.004 


-0.2 


<1.4 


HD /57U 


4.4±0.5 


0.004 


0.3 


<9.7 


1 n 7^n n 

iu. ( ±u.y 


0.082 


4.3 


0.8 


HU 9ozb 


-7.4±0.9 


-0.003 


-0.3 


<9.7 


-1.3±1.5 


-0.004 


-0.2 


<0.5 


TT n -i nonn 
H4J lUoUU 


8.8±0.3 


0.015 


1.2 


<13 


-1.6±0.8 


-0.024 


-1.1 


<0.7 




T 7 A _i_n A 


0.010 


0.8 


<13 


0.0±3.3 


-0.006 


-0.1 


<0.7 


tiu 4ooo4 


0.1±0.6 


-0.001 


0.0 


<13.9 


0.5±0.9 


0.004 


0.2 


<0.7 


tiu 555/5 


4.1±0.5 


0.006 


n a 
U.4 


<11.1 


0.8±0.9 


0.008 


n a 
U.4 


<0.6 


HJJ 5oo55 


3.1±0.5 


0.005 


n a 
U.4 


<9.7 


-0.4±0.9 


-0.004 


-0.2 


<0.5 


TJn fiOQQn3 

44 4J oyooU 


72.3±4.5 


0.099 


T n 


41.6 


•4z±4 


0.50 


24.4 


8.0 


un wnnc; 


OO Q_L1 O 

ZZ.o±4.z 


0.028 


2.2 


9.5 


7.0±0.5 


0.083 


4.5 


1.1 


HU 757o2 


1 A a _i_n c 
44.4±U.5 


0.019 


1.5 


<13.9 


1.5±0.9 


0.019 


0.9 


<0.7 


TT n T'C 1 CI 


i n n_i_n q 
4U.U±U.o 


0.018 


1.4 


<13 


10.1i0.7 


0.151 


( .6 


2 


xlJJ o4 / o f 


oi 1 _i_n 
Z4. 4±U.D 


n ni i 
U.U4 ( 


1 A 
4.4 


<44.4 


n 1 1 
-U.v5±4. 4 


n nno 
-U.UUz 


n i 
-U. 4 


<U.D 


riJJ ob ( zo 


-4o. i ±U. / 


n ni o 
-U.U4Z 


n n 

-u.y 


<4o 


n i i n q 
U. 4±U.o 


n nno 
U.UUz 


n i 
U.4 


<U. ( 


un ac; 1 no 
tiu y5lzo 


oo 1 _i_n n 

-zo.i±u.y 


-0.022 


1 T 

-1.7 


<11.1 


O O 1 1 Q 


n nnc 
U.UUb 


0.3 


<U.D 


un i m trm 4 
41 JJ 1U45U1 


-yo.4±4.o 


n 071 
-U.U f 4 


-5.D 


^ 1 Q O 


U.D±4.4 


n nnt; 
U.UU5 


n q 
U.o 


<U. / 


u n 1 om 9c 
H4J IzUloD 


45.o±U.b 


U.U44 


n q 
U.o 


<y. / 


1 o i 1 a 


n nnv 
-U.UU f 


n q 
-U.o 


<U.O 


un 1 oo q 1 1 
H4J Izooll 


-0.9±0.4 


-0.003 


-0.2 


<16 


1.4±0.6 


0.037 


1.6 


<0.8 


u n 1 oonno 
H4J looUUz 


1.8±0.6 


0.002 


0.1 


<11.1 


2.9±1.0 


0.025 


1.3 


<0.6 


u n 1 Q£tr\iz a 
tiu 4ooUo4 


-6.6±0.6 


-0.007 


-0.6 


<9.7 


0.2±1.0 


0.002 


0.1 


<0.5 


Un 1 A OQ7Q 
444J 44ZO 1 O 


1 o_Ln i 


U.UU4 


n n 
U.U 


<44.4 




n m v 
-U.U4 f 


n q 
-U.o 


<U.D 


44 iJ 44ozoo 


-28.1il.2 


-0.032 


-2.5 


<13 


-1.6i0.9 


-0.018 


-0.9 


<0.7 


u n 1 C/inoo 
414J 454Uoo 


-5.3±0.3 


-0.013 


-1.0 


<13.9 


-0.6il.l 


-0.017 


-0.6 


<0.7 


44 iJ 4bObZU 


-0.3±0.4 


0.000 


0.0 


<26.8 


-0.2i0.6 


-0.002 


-0.1 


<1.4 


44 IJ 4D5454 


45. i ±U.O 


U.UlO 


1.3 


<8.2 


-0.9i0.8 


-0.009 


n c 
-U.5 


<U.4 


HD 173667 


46.4±0.9 


0.023 


1.8 


<9.7 


9.7±1.5 


0.043 


2.3 


<0.5 


HD 181321 


0.3±0.4 


0.001 


0.1 


<13 


0.0±0.6 


0.002 


0.1 


<0.7 


HD 185144 


-5.5±2.1 


-0.003 


-0.2 


<16 


3.2±2.0 


0.011 


0.6 


<0.8 


HD 188376 


12.8±1.1 


0.005 


0.4 


<13 


-3.5±2.2 


-0.013 


-0.7 


<0.7 


HD 191408 


-8.7±1.0 


-0.004 


-0.3 


<26.8 


4.4±1.8 


0.020 


1.0 


<1.4 


HD 196378 


-15.3±1.2 


-0.012 


-0.9 


<9.7 


2.1±0.8 


0.016 


0.8 


<0.5 


HD 203608 


-38.2±0.8 


-0.016 


-1.3 


<9.7 


-5.2±1.4 


-0.019 


-1.0 


<0.5 


HD 206860 


-0.6±0.4 


-0.001 


-0.1 


<11.1 


5.3±0.7 


0.085 


4.0 


0.9 


HD 212330 


16.7±0.8 
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-0.2 


<0.7 
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0.003 
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<13 
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0.002 


0.2 
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-0.5 
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14±6 








0.74±0.31 



1 Limits on Ld us t/L* from using Eqn (2) and three times the dispersion in the fractional excess of the whole population. 
2 Excess divided by dispersion averaged over entire sample. 
3 Data on excess taken from Beichman ct al 2005b. 
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4 Star shows significant flat field errors after re-calibration and is not considered further. 
5 Average value of Ldust/^* (3a limits or detections) after 2a rejection of outliers. 
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Table 3. New MIPS Observations 1 



Star 


AOR 


F v (24 /nil) 2 
(mjy) 


Phot. 3 
(mjy) 


Obs./Phot. 


F u (70 fxm) 
(mjy) 


Phot, 
(mjy) 


Obs./Phot. 


X70 4 


Ld/L-t (70 fim) 


HD693 


4030976 


254 


265 


0.96 


37.7±5.6 


30.0 


1.26 


1.37 


< 0.4 x 10" 5 


HD9826 5 


4033280 


519 


539 


0.96 


51.2±10.3 


60.7 


0.84 


-0.92 


< 0.4 x 10" 5 


HD154088 


4061696 


87 


70 


1.24 


-0.1±6.6 


7.8 




-1.19 


< 2.7 x 10" 5 


HD191408 


4055040 


443 


435 


1.02 


46.7±12.2 


50.1 


0.93 


-0.28 


< 1.5 x 10" 5 


HD196378 


4055808 


227 


254 


0.90 


30.5±5.6 


28.9 


1.06 


0.29 


< 0.4 x 10~ 5 


HD212330 


4058368 


235 


234 


1.01 


26.4±5.8 


26.5 


1.00 


-0.01 


< 0.7 x 10" B 


HD217014 


4058880 


188 


178 


1.06 


25.2±5.3 


20.0 


1.26 


1.00 


< 0.2 x 10~ 5 


HD217813 


4062464 


59 


70 


0.85 


3.7±4.8 


7.9 


0.47 


-0.87 


< 1.8 x 10~ 5 


HD222368 


4060160 


514 


534 


0.96 


54.±10.2 


60.4 


0.91 


-0.56 


< 0.4 x 10" 5 


Average 








1.00±0.04 






0.92±0.09 







1 MIPS data obtained after data presented in Beichman et al 2005a and Brydcn 2006 
Statistical noise is negligible. 

3 Photospheric prediction based on extrapolation from 2MASS and shorter wavelength data as described in the text and in 
Bryden et al. 2006. 

4 X70 defined as (Observed — Photosphere) /Noise. 

5 Ups And 



-25 - 



Table 4. MIPS Excesses 



Star F u (70 /im) Excess (mJy) 1 Ldust/-^* (70 /xm) Ld us t/L± (10-70 /im) 

HD 7570 14±6 0.38xl0~ 5 0.43xl0~ 5 

HD 69830 2 <3 < 0.5 x 10~ 5 20xl0~ 5 

HD 76151 16±4 l.OxlO" 5 1.2xl0~ 5 

HD 72905 3 24±4 1.2xl0~ 5 2.9xl0~ 5 

HD 128311 18±4 2.6xl0~ 5 

HD 206860 15±4 0.92xl0~ 5 15xlQ- 5 



^ryden et al 2006 
2 Beichman et al 2005b 

3 The star HD72905 has a 70/zm flux density of 24±3 mJy (Bryden et al. 2006), a factor of 
2-3 lower than measured by ISO-PHOT at 60 and 90 /im (Spangler et al. 2001). We adopt 
the MIPS value. 
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Table 5. Model Fits to Spectra 



Star 


Rbb 
(AU) 


(K) 


R1-R2 
(AU) 


Ti-T 2 

(K) 


Optical Depth 
t (10 /im) 


Reduced 

x 2 


M e 


M e <10 km 


HD 7570 
HD 76151 
HD 128311 
HD 206860 


7.8±0.8 
8.1±0.6 

> 15 
13.3±1.2 


98±20 
90±14 
< 55 
73±13 


11-12 
7.6-8.7 

13.2-14.3 


83-86 
93-86 

74-71 


1.4xl0~ 5 
21xl0~ 5 

21xl0 -5 


0.68 1 
0.62 1 

0.5 1 


1.3xl0" 6 
1.6xl0~ 6 

2.7xl0~ 6 


0.042 
0.052 

0.087 


HD 72905 

0.25 (im silicate 
and 




0.03-0.43 


1500-540 


10xl0- 5 




7.2xl0- 1() 


1.4X10- 4 


0.25 (im silicate 

or 

10 fim silicate 






109-315 


55-36 


33 xl0~ 6 


0.83 2 


HOxlO" 6 


22 


14±1 


70±10 


12.2-15.9 


67-63 


7xl0~ 5 


0.83 2 


3.3xl0~ 6 


0.10 



1 X > 21/tm, 74 degrees of freedom 
2 A > 8/xm, 256 d.o.f. 
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Fig. 1. — Histograms showing the distribution of fractional excess ((Observed- 
Photosphere) /Photosphere) estimated within pass bands at 8.5-13 (left) and 30-34 /im (right). 
The point for HD 69830 is off-scale to the right at 30-34 /iin with a relative excess of 0.5 
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Fig. 2. — These plots show the IRS spectrum (top), excess flux (middle), and fractional ex- 
cess (bottom) as a function of wavelength for HD128311, HD 216437, HD 206860, HD7570, 
HD 72905 and HD76151. The first two stars have no measurable excess and are shown as 
representative of the majority of the sample. Also plotted in the lower figure are the approx- 
imate 2a dispersion in the deviations from the photospheric models. Fractional deviations 
between these lines should be regarded with skepticism. 
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Fig. 2. — continued. 
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Fig. 2. — continued. 
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Fig. 2. — continued. 
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Fig. 2. — continued. 
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Fig. 2. — continued. 
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Fig. 3. — Histograms showing the detections or (mostly) 3<r upper limits on L^ ust /L^ for dust 
emitting in 30-34 /zm band (dashed lines) and the 8.5-13 /zm band (solid lines) for the stars 
in this survey. These limits take into account the effective temperature, T e jj, and fractional 
excess for each star as described in the text and are given in units of the solar system's 
approximate level of zodiacal emission, L dust /L* = 1CT 7 (Backman & Paresce 1993). The 
positions of the positive detections are indicated along with the star's name. An —8 or —30 
is appended to the name to resolve any ambiguity as to which band is being shown. The 
remainder of the values are upper limits. 
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Fig. 4. — IRS spectra of those stars with excesses (solid line) along with their best fitting 
models (dashed line). The offset in SL-1 spectrum of HD 76151 is not statistically significant 
in terms of its fractional excess(1.4cx). The spectrum of HD 69830, with its prominent IRS 
excess, is presented in Beichman et al. 2005b. 
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Fig. 5. — For each star with planets the minimum and maximum extents of the gravitational 
influence of that star's planets are plotted. The lack of emission in the IRS wavelength region 
is consistent with a dearth of material interior to about 5 AU. A partial explanation may be 
the influence of planets on dust in disks. 



